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ABSTRACT 

We report on the analysis of high-precision space-based photometry of the roAp (rapidly oscillating Ap) stars HD9289, HD 99563, 
and HD 134214. All three stars were observed by the MOST satellite for more than 25 days, allowing unprecedented views of their 
pulsation. We find previously unknown candidate frequencies in all three stars. We establish the rotation period of HD 9289 (8.5 d) 
for the first time and show that the star is pulsating in two modes that show different mode geometries. We present a detailed analysis 
of HD 99563 's mode multiplet and find a new candidate frequency which appears independent of the previously known mode. Finally, 
we report on 1 1 detected pulsation frequencies in HD 134214, 9 of which were never before detected in photometry, and 3 of which 
are completely new detections. Thanks to the unprecedentedly small frequency uncertainties, the p-mode spectrum of HD 134214 can 
be seen to have a well defined large frequency spacing similar to the well-studied roAp star HD 24712 (HR 1217). 

Key words, stars: oscillations ~ stars: chemically peculiar - stars: magnetic fields ~ stars: rotation 



1. Introduction 

Rapidly oscillating Ap (roAp) stars are cool magnetic chemi- 
cally peculiar stars with spectral types from early A to FO, lumi- 
nosity class V. They were discovered by Kurtz | ( 1982| l and pul- 
sate with low amplitudes (< 13mmag) at periods between 5 to 
21 minutes. The oscillations are due to high-overtone (« > 15) 
non-radial p-mode pulsations. The oblique pulsator model (here- 



after OPM) as described by Kurtz ( 1982 1 has been, to a large ex 



tent, a very applicable explanation of their observed amplitude 
variability. It predicts that the pulsation axis is aligned with the 
magnetic axis, which itself is oblique to the rotation axis of the 
star. The OPM has since been refined to include additional ef- 
fects like the Coriolis force ( Bigot & Dziembowski|2002) l. For 
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revie ws on primarily observatio nal aspects of roAp star s we re- 
fer to [Kurtz & Martinez l(t2000l) and IKochukhov' (20071), and on 
theoretical aspects, to|Shibahashi (2003 ) and Cunha (2005). 

The roAp stars offer great potential for asteroseismology. 
This potential is difficult to realise due to the small number of 
roAp stars known, their short pulsation periods and low am- 
plitudes, and (apparently) very few observable modes. Recent 
investigations with the WIRE, MOST, and Kepler satellites, as 
well as groundbased high-resolution spectroscopy, have allowed 
us to study these objects with an unprecedented level of detail. 

Encouraged by the success of previous MOST observations 
of roAp stars y Equ ( [Gruberbauer et al.|2008| l and 10 Aql ( [Huber[ 
[et al.|20(r8( l, new observations of other roAp stars accessible to 
MOST were planned and peiformed with the primary goal of de- 
tecting additional frequencies. In this paper we report our find- 
ings on the targets HD 9289, HD 99563, and HD 134214. 

HD 9289 (V ^ 9.4) was orig inally classified as ApScEu by 
Bid elman & MacConnell ( 1973|l and is now considered to be an 
A3p star (see Ochsenbein et al.[[1981[ and references therein). 
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Table 1. Summary of the MOST data sets of HD 9289, HD 
99563, and HD 134214. 



HD 




Mo, 


'red 


A'red 


Ccff 




[d] 




[d] 






HD 9289 


26.27 


41884 


26.06 


37635 


0.632 


HD 99563 


27.53 


62261 


26.68 


55176 


0.773 


HD 134214 


25.96 


35099 


23.92 


31244 


0.459 



Notes. /toc> Mot^ total timebase and number of data points for MOST 
observing run; f,ed, Med: timebase and number of data points for final 
light curve; c^ff: elfective duty cycle for final light curve (see text for a 
definition) 



It was discovered to be a rapidly oscillating Ap star by Kurtz 
et al. ( 1994) 1. The authors observed HD 9289 on 22 nights, span- 
ning several months from 1993 to 1994, and identified three in- 
dividual frequencies around ~ 1550 yuHz. The latest results on 
this star were obtained spectroscopically by Ryabchikova et al. 
( |2007| l, who were able to refine the values of some of its fun- 
damental parameters. Most notably, they were able to estimate 
vsin/= 10.5 km s"' and (Bp) = 2.0 kG for the first time. 

The first discovery of roAp pulsation in HD 99563 was re- 
ported by Dorokhova & Dorokhov (1998i. The star was later 
studied extensively using long-timebase, ground-based photom- 
etry (Handler et a l. 2006), as well as spectroscopy ( Fi'eyh ammer| 



et al. 2009| l. HD99563 is particularly interesting for the latter 



type of study because its pulsation shows large radial velocity 
amplitudes compared to other objects of this class. In all stud- 
ies conducted so far it was found to pulsate in a single intrinsic 
pulsation mode which is split into a multiplet due to the effects 
predicted by the OPM. 

HD 134214 was initially recognised to be a roAp star in 
1985. A first thorough frequency analysis resulted in a single 
prominent mode with an amplitude of several mmag, at a fre- 
quency of 2949.6 (iHz ( Kreidl & Kurtz|198 6). No signs of rota- 
tional multiplets were detected. Subsequent observations did not 
reveal additional frequencies, but indicated that the frequency of 
the mode was periodically changing with a period of about 250 
days ( |Kreidl et al.|1 994). More than a decade later, time-resolved 
high-resolution spectroscopy of several roAp stars revealed sig- 
natures of previously unseen frequencies in almost all the stars 
in the sample, including HD 134214 ( [Kurtz et alJ|[2006] ). This 
new variability was only observed in specific lines in the upper 
stellar atmosphere. Apparently, their amplitudes were too small 
for them to be detected with previous ground-based photometric 
observations. A 10-hour run with the MOST satellite could not 
confirm any of the newly reported signals, but it set an upper 2-cr 
detection limit of 0.36 mmag for the corresponding amplitudes 
in the MOST bandpass ( Cameron et al.|2006|l. A later spectro- 
scopic run dedicated to HD 134214 ( Kurtz et al.|2007| l, was in- 
tended to clarify the values of the newly seen frequencies. Due to 
the relatively short duration of the observations, and the apparent 
dependence of the frequency solution on the used spectral lines, 
no general list of new modes could be given. Nonetheless the 
general properties of the previously detected signals was con- 
firmed and one of the additional frequencies detected in photo- 
metric observations. 



2. Observations and data reduction 

MOST (Microvariability and Oscillation of STars) is a Canadian 
microsatellite designed for the rapid-cadence observation of stel- 
lar photometric variability on a level of a few parts-per-million 



(ppm). The satellite is in a sun-synchronous polar orbit, above 
the terminator, and has an orbital period of 101 min. It is able 
to observe bright stars in its continuous viewing zone (CVZ) for 
up to 2 months. It houses a Maksutov telescope with an aper- 
ture of 15 cm, a custom broadband filter (3500 - 7000 A), and 
two CCDs. One CCD was originally designed to be exclusively 
employed for satellite guiding (Tracking CCD), while the other 
CCD should collect the scientific observations (Science CCD). 
To accomplish its task more accurately, the latter is equipped 
with a Fabry -lens array covering a small part of the detector area. 
When observing stars in the upper brightness range of the instru- 
ment, the Fabry lenses are used to reduce the jitter of the image 
that is caused by pointing instabilities and to avoid saturation by 
spreading the light over a very large number of pixels. In con- 
trast, dimmer programme stars and guide stars are observed in 
the open field (Direct Imaging - DI). A more detailed description 
of the early mission is given by Walker et al. (2003 ). In 2006 the 
Tracking CCD stopped working, and the Science-CCD software 
was modified to additionally provide the information for the at- 
titude control system (ACS). Since the ACS constantly needs 
information about the satellite pointing, stars are observed with 
short exposure times, but the individual exposures are stacked to 
reach and improve upon the original levels of photometric pre- 
cision achieved by MOST. 

As explained by |Reegen e t al.'('2006), MOST suffers varying 
amounts of stray light as it orbits about the partially illuminated 
Earth, which causes artificial signal at multiples of the MOST or- 
bital frequency. In addition, this non-stellar signal is modulated 
throughout an observational run because of changes in the aver- 
age albedo of the illuminated limb of the Earth as seen from the 
position of the satellite. This modulation occurs on timescales 
comparable to a sidereal day. Longer-term instrumental drifts 
are sometimes also present, but MOST has shown instrumen- 
tal stability of about 1 mmag over intervals of years. The MOST 
team has become very familiar with these effects, and photomet- 
ric reduction pipelines for the different MOST observing modes 



(|Reegen et al ||2006[ [RoweeraL][2006l [Huber & Reegen||2008 
Har eter et al.|2008| effectively filter the artifacts in many cases 
while leaving the stellar signals intact. 

All three roAp stars discussed in this paper have been ob- 
served using the DI mode. The general properties of these data 
sets are summarised in Table[T] This table also includes what we 
call the "effective duty cycle" 



Ceif 



m(N-i) 



tn 



(1) 



where (Af) is the median time stamp difference between two ad- 
jacent data points, is the total number of observations, and 
fobs = t (N) - t (0) is the total time baseline of the resulting time 
series. The median is not very sensitive to outliers, and so closely 
reflects the effective sampling of the time series, even if sporadic 
under-sampling occurs due to larger gaps. 

One can also define an effective sampUng frequency 

similar to the definition of the Nyquist frequency. This, in combi- 
nation with the window function, is a good measure of the actual 
sampling quality of the data. 

2.1. HD9289 

MOST observed HD 9289 for about 22 of the 26 days during 
September 27 to October 23, 2007. After the first day of obser- 
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Fig.l. Reduced MOST light curves of (a) HD9289, (b) HD 99563, and (c) HD 134214. The inserts shows a small portion of a 
version of the light curve where any variability not due to roAp pulsation has been subtracted. For the insert in panel (a), a 5-point 
running average is displayed instead of the unbinned data, in order to enhance the clarity of the weaker pulsation signature. Note 
that JD2000 = 2451545.0 



vation a 4-day gap occurred, and during the first 16.4 days of ob- 
servation MOST alternated between two fields of view. During 
this setup HD 9289 was monitored for 4 consecutive orbits, after 
which MOST switched to the other primary target for the dura- 
tion of 1 orbit, and then back to HD 9289. For the remaining 10.6 
days, however, HD 9289 was observed without interruption. The 
vast majority (> 99 %) of the measurements consisted of image 
stacks of 20 exposures each, with a total exposure time of 32.13 
seconds. 

The data were reduced with the Dl-pipeline developed by 
Huber & Reegen (2008). Different pipeline parameters were set 



to produce two versions of the light curve: the "final light curve" 
version, where instrumental effects are suppressed as far as pos- 
sible, and the "consistency check" version that underwent a less 
invasive treatment and more closely reflects the raw data. For 
each version of the reduction, about 4200 frames (about 10% of 



the total number in the time series) were rejected during the ini- 
tial stages of the reduction process because of heavy stray light 
pollution, and/or excessive cosmic -ray incidence in the proxim- 
ity of the Southern Atlantic Anomaly (SAA). The values cited 
in Table[T]refer to the light curve with the smallest instrumental 
contributions. Apart from increased scatter at high-stray light or- 
bit phases and possible long-term trends, no obvious instrumen- 
tal effects are present in the final light curve (see Fig.[T]i. It shows 
a regular double-hump variation with a period of about 8.5 days 
which has been tracked for almost 3 cycles. This is actually rota- 
tional modulation of the mean stellar brightness, probably due to 
spots. It is the first time that such a signature has been detected 
for HD 9289. The insert in panel (a) of Fig. [T] also demonstrates 
the presence of the short-period roAp pulsations. 
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Fig. 2. Light curve of HD 99563 after processing with the reduc- 
tion pipeline (top) and after subtraction the instrumental long- 
term trends (bottom). The dashed vertical lines separate the in- 
dividually fitted subsets of the time series. 



Fig. 3. Same as Fig.jijbut for HD 134214. On top of the long- 
subset polynomial fits the top panel also shows a 200 pt running 
average used as a high-pass filter. 



2.2. HD 99563 

MOST monitored HD 99563 for about 26 days during February 
25 to March 23, 2008. However, the first day of observation did 
not yield any usable data. A second target was observed in a 
different field, hence MOST was re-oriented every 3rd orbit for 
about 50 minutes, decreasing the duty cycle. Almost all of the 
remaining frames consist of 15 stacked exposures, with a total 
integration time of 32.32 seconds. With respect to the reduction, 
the same steps were performed as for HD 9289. However, due 
to a change in the observation setup at JD - JD2000 ~ 2992.5 
(JD2000 = 2451545.0) there is a half-day gap in the time series. 
The data obtained beyond this gap are of slightly lower quality 
due to an increase in stray light, and a long-term instrumental 
trend, which is also reflected in the board temperature teleme- 
tr^ To improve the situation the latter parts of the light curve 
were treated separately using different parameters for the reduc- 
tion pipeline. The data set was then homogenised by fitting and 
subtracting polynomials of second order to the two parts (see 
Fig.|2]i. The final light curve is shown in Fig.[T]3. There is obvi- 
ous variability on timescales much longer than what is expected 
from the pulsation alone. Some of the variability is caused by the 
OPM effect, which modulates the pulsation amplitude over one 
stellar rotation cycle. A second variability component is due to 
rotational modulation of the mean stellar brightness, which has 
already been detected in previous studies (e.g., [Handler et al.| 
2006| . Because of the broad MOST bandpass (about 350 - 700 
nm) and the wavelength dependence of Ap rotational modula- 
tion, this signal is much weaker than what Handler et al. found 
in the ground-based data with a standard blue filter It is also less 
pronounced than the spot effects in HD 9289. A third component 
is purely instrumental and caused by complex onboard tempera- 
ture variations due to the target switching. We did not attempt to 
remove this component during the reduction process in order to 
not distort the intrinsic rotational effects on the pulsation. 



' The CCD focal plane temperature is kept constant to within 0. 1 K, 
but the rest of the satellite temperature can vary by about 2 K during 
each MOST orbit. This effect can be more pronounced if MOST is 
switching between two different, widely separated fields during each 
orbit 



2.3. HD 134214 

HD 134214 was observed for about 26 days during April 18 to 
May 14 2008. Each MOST orbit was shared equally between this 
star and another science target field throughout the run, reduc- 
ing the duty cycle to about 50 per cent. The majority of images 
are stacks of 19 exposures (total integration time: 30.2 seconds). 
Again, we used the same steps to reduce the data as in the anal- 
ysis of HD 9289, yielding two versions of the light curve. The 
only difference is that, after reduction, long-term trends, which 
correlate with board temperature variation, are still visible in the 
data. These have been removed from the light curve (see Fig.[3]l 
in two ways. First we fitted and subtracted fourth order polyno- 
mials from three distinct subsets of the light curve clearly sep- 
arated by larger gaps. However, we also noticed a modulation 
of the board temperature on short time scales due to the target 
switching. Due to the absence of any obvious rotational signa- 
ture, we therefore also used a high-pass filter, removing any vari- 
ations on time-scales longer than the MOST orbital period. The 
latter approach produced fewer instrumental artefacts, and the 
result is referred to as the "final light curve" for HD 134214. 
However, we repeated all further analysis steps also for the less 
intrusive approach to make sure that we did not actually intro- 
duce any fake signal. The final light curve is dominated by the 
roAp pulsation and does not show any significant long-term vari- 
ability after removing the instrumental trends. 



3. Frequency analysis 

The final light curves were analysed using the latest version of 
SigSpec ^Reegen 2007 1. SigSpec 2.0 performs an iterative multi- 
cosine fit to time-series data, based on the "spectral significance" 
(or sig) as a significance criterion. The sig is defined as the neg- 
ative logarithm of the probability that a periodic signal of the 
form a cos (27rvf - 0) (where a is the amplitude, v is the fre- 
quency, and 6 is usually called the phase) is caused by white 
noise. SigSpec starts by computing the sig for a given range in 
frequency and finding the frequency, amplitude, and phase that 
correspond to the smallest false-alarm probability. It then refines 
the solution by performing a least-squares fit in the time domain. 
Subsequently, the solution is subtracted from the time-series, and 
the sig is computed for the residuals. This continues iteratively. 
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each time adding another component to the least squares multi- 
cosine fit, until no signal is found with a sig above a certain 
threshold. Formally, an amplitude-to-noise ratio of 4, which is 
generally assumed to be a conservative lower limit, corresponds 
to sig ^ 5.46. Self-evidently, the definition of the sig is only 
truly valid in the case of a single sinusoidal signal in the pres- 
ence of white noise. The results therefore have to be interpreted 
with care, and often a higher threshold than 5.46 is sensible. 

Even when data reduction is performed thoroughly, a number 
of significant instrumental signals due to stray light or target field 
switches can remain in MOST time series. Any gaps due to target 
switching and/or outlier removal during high stray light phases 
make it imperative to correctly pre-whiten these signals together 
with the intrinsic frequencies. To properly take this into account 
in the analysis of the three roAp stars we proceeded as follows. 

In the first step, SigSpec was used to calculate the most prob- 
able multi-cosine solution down to a sig threshold of 5.46 for 
each reduced light curve. The upper limit in frequency for the 
SigSpec analysis was set for each star so we could detect a pos- 
sible harmonic to the main pulsation frequency. The same cal- 
culation was performed for time series composed of the mean 
background measurements, and also board temperature mea- 
surements, as provided by the reduction pipeline. For each star, 
the resulting frequencies of all time series were then compared 
and put into one of three categories. 

1. Instrumental artefacts, (a) Identifiable in both the stellar 
photometry and in either background or temperature teleme- 
try, and/or (b) obviously related to a harmonic of MOST or- 
bital frequency 14.19 d"'). Their amplitudes vary drasti- 
cally in the stellar time series when different pipeline param- 
eters are used. 

2. Intrinsic low-frequency variability. Any highly significant 
signals at lower frequencies (i.e., due to rotational modula- 
tion caused by spots) and found in both reduced light curves 
but not in the background/telemetry. 

3. Intrinsic pulsation (candidates). Possible pulsation fre- 
quencies, in particular signals in the range associated with 
roAp oscillations, that are found in the more and the 
less invasively reduced light curves but not in the back- 
ground/telemetry. 

Frequencies in categories 2 and 3 were checked for consis- 
tency with Period04 (Lenz & Breger [200 5). For further time- 
resolved analysis of the pulsation candidates, all signals except 
those in category 3 were subtracted from the light curves. The 
uncertainties for frequencies, amplitudes, and phases are ana- 
lytical lower limits calculated with Period04. We cross-checked 
these values with the SigSpec results using the heuristic rela- 
tions from [Kallinger et aT] ( |200 8 ). The latter provide a reliable 
upper limit in the uncertainties. Interpreting these as 4cr confi- 
dence limits we find both results to be highly consistent. 



3.1. HD9289 

The results of the frequency analysis for HD 9289 are presented 
in Table |2] The MOST data reveal the presence of 6 intrin- 
sic significant frequencies between 1570 and 1590 ^iHz. In the 
low-frequency range, we can also find signals corresponding to 
the rotational modulation. However, in contrast to the high fre- 
quency signals, the exact values of these low frequencies depend 
on the reduction procedure used. This reflects the fact that it is 
very difficult to preserve the exact nature of long-term trends 
in space photometry due to the lack of calibration sources and 
procedures. This can lead to problems with analyses using the 



Table 2. Frequency analysis results for HD 9289. 



id 


frequency 


amplitude 


e 


sig 


SNR 




[nHz] 


[mmag] 


[rad] 






Vl 


1584.936(7) 


0.63(2) 


5.05(2) 


259.8 


28.6 


Vi 


1582.22(2) 


0.24(2) 


4.68(6) 


48.4 


10.8 


V3 


1587.67(2) 


0.16(2) 


5.8(1) 


20.4 


7.1 


V4 


1573.27(3) 


0.15(2) 


6.0(1) 


23.2 


6.5 


vs 


1571.92(3) 


0.14(2) 


1.6(1) 


16.4 


6.1 


Vb 


1577.33(4) 


0.11(2) 


5.7(1) 


10.5 


4.9 


^■rot.PDM 


1.35(1) 










^rot.puls 


1.361(7) 











Notes. Phase is given according to the definition cos {2nvt - 0) with 
respect to the epoch HJD 2454371.05691. The numbers in brackets are 
the errors in units of the last digit. SNR stands for signal-to-noise ratio. 



Discrete Fourier transform (DFT), which expects a periodic har- 
monic behaviour of the variations that are analysed. 

Since we can assume that the rotational modulation of roAp 
stars, caused by spots, is stable across many cycles, any devia- 
tions are likely instrumental. We can therefore apply the phase 
dispersion minimization technique (PDM, Stellingwerf 1 1 978 1 to 
obtain an estimate of the rotation frequency based on the period- 
icity of the signal, more independent of the exact preservation of 
its shape. We used the PDM technique on the two reduced light 
curve versions, as well as on a raw data version that was filtered 
as to best preserve the intrinsic long term variations. We then 
calculated the mean and standard deviation of the resulting rota- 
tion frequencies, noted as Vrot.poM in Table|2] Luckily, this crude 
estimate can be improved upon, since there are other ways to get 
a less reduction-dependent estimate of the rotation frequency, 
^rot.puis, as discussed in Sect. 



4.1 



A large number of additional weak signals with frequencies 
below 10 d"' are present as well, but a comparison with the back- 
ground and board temperature reveals that most of them are in- 
strumental. We conclude that we cannot reliably assess the in- 
trinsic nature of any additional low-frequency variability. Fig.|4] 
compares the amplitude spectrum before and after pre- whitening 
all signals that are members of category 1 and 2. Fig. |7] shows a 
schematic spectrum of the pulsation multiplets. 



3.2. HD 99563 

As presented in Fig. [5] a large number of instrumental arte- 
facts with substantial amplitudes can be identified in the low- 
frequency regime. This stems from target switching, as dis- 
cussed in Sect. 2.2 Nonetheless, we can identify what seems 
to be the stellar rotation frequency and one of its harmonics. 
However, the value found for the "fundamental" is 3.902(5) [xHz, 



which is inconsistent with the precise value given by 'Handler 
et al. (2006|. The second frequency interpreted to be the first 
harmonic found at 7.970(3) ^iHz is in better agreement with the 
previously published value, but it is still too far off given the for- 
mal uncertainties, and incompatible with our value of the "fun- 
damental". A PDM analysis of the least noisy parts of the data 
yields a value of 3.95 ^iHz, but the result is not well constrained 
due to the small amplitude of the variation. We thus conclude 
that we cannot improve the accuracy of the rotation period using 
the mean light variation. 

Concerning pulsation frequency candidates, we treated 
HD 99563 in a special way, since its main (and hitherto only) 
pulsation frequency is already so well studied. We first per- 
formed our unbiased, unconditioned analysis as explained at the 
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Fig. 4. DFT of the HD9289 data without any pre-whitening 
(grey) and after pre-whitening of all category 1 and 2 signals 
as defined in Sect.|3](black). The insert shows the window func- 
tion of the data set. Vertical dotted lines show the MOST orbital 
frequency and its integer mulitples. 



beginning of the section. The results for this approach are pre- 
sented in Table|3] As expected, we find 7 more or less equally 
spaced frequencies centred at » 1557.7 ^iHz. We also find a sec- 
ond group of frequencies composed of integer multiples (com- 
monly called "harmonics") and linear combinations of the sep- 
tuplet frequencies. The broad picture is mostly consistent with 
previous observations. We do detect the (-2,H-2) components 
that Handler et al. could not find in their ground-based data. 



However, we cannot find the outermost (-4,H-4) nonuplet compo- 
nents that|Fre yhammer et al.| report. A closer look reveals addi- 
tional discrepancies. The frequency value for the central compo- 
nent of the septuplet differs by almost 3cr. We also find a number 
of additional significant frequencies in between the septuplet. 
These candidates are labeled vg to vn in Table[3] 

Most of the remaining frequencies found are consistent with 
well-known stray light artefacts, situated at multiples of the 
MOST orbital frequency. What is left can be explained through 
another effect: modulation of the stellar signal with the MOST 
orbital frequency vmost- We find formally significant frequen- 
cies at V2 + «vmost and V3 ± «vmost, where n is an integer up to 
« = 3. This is not aliasing, as one might think when looking at 
the window function of the data set, but rather an actual instru- 
mental modulation of the measured stellar signal. These effects 
cannot be easily corrected for, but they can easily be identified 
in the course of the frequency analysis. 

The rotation period of HD 99563 can't be determined to high 
accuracy from the light modulation signal in the MOST time se- 
ries. But it can be obtained independently, and in an unbiased 
way, from the average spacing of the 0PM multiplet. We find 
(vrot) - 3.985(7) [xHz, which is compatible with the literature, 
but obviously not precise enough to be considered an improve- 
ment. 

We also chose a second approach that would more closely 
reflect our prior knowledge of HD 99563's pulsation and rota- 
tion. Handler et al.j give a quite precise value for the rotation 
frequency of the star v,ot,Handier - 3.9749(1) [.iHz. Using this 
value, we forced a fit to a nonuplet with the corresponding spac- 
ing, looking for whatever frequency gave the best results. The 
fitting procedure was carried out using Period04, after all cat- 
egory 1 and 2 signals had already been subtracted. The result 
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Fig. 5. Same as Fig. [4] but for HD 99563. Notice the presence 
of a small signal "hump" due to the harmonic to the main pulsa- 
tion frequency, and combination frequencies, at about 3110 [xHz. 
Another small hump at about 2380 (-iHz is caused by aliasing. 



Table 3. Unconditional multi -frequency fit results for HD 99563. 



id 




frequency 


amplitude 


e 


sig 


SNR 






[nHz] 


[mmag] 


[rad] 






Vi = 




1557.68(1) 


0.23(1) 


4.58(4) 


88.1 


15.8 


V2 ~ 


Vp-l 


1553.683(2) 


1.50(1) 


2.634(7) 


1969.5 


102.5 


V3 a; 


Vp+l 


1561.657(3) 


1.11(1) 


4.796(9) 


1184.8 


73.3 


V4 a; 


Vp-2 


1549.68(4) 


0.06(1) 


2.8(2) 


6.7 


4.2 


V5 a; 


Vp+2 


1565.67(3) 


0.08(1) 


0.6(1) 


11.9 


5.3 


Vf, a 


Vp-3 


1545.76(2) 


0.16(1) 


1.08(6) 


44.0 


10.8 


v^ a; 


Vp+3 


1569.56(3) 


0.10(1) 


0.4(1) 


18.5 


6.5 


Vg 


1553.21(3) 


0.12(1) 


2.5(1) 


10.9 


8.3 


vg 




1561.83(4) 


0.10(1) 


1.0(1) 


7.7 


6.9 






1552.35(4) 


0.06(1) 


2.5(2) 


8.3 


4.2 


Vil 




1554.01(4) 


0.17(2) 


5.42(5) 


6.4 


12.5 


V12 = 


= I'cand 


1555.84(4) 


0.06(1) 


0.6(2) 


7.1 


4.2 


2vp 




3115.36(2) 


0.13(1) 


2.78(8) 


32.7 


9.7 


2yp_ 




3107.46(3) 


0.09(1) 


0.3(1) 


18.1 


7.1 


2vp - 




3119.38(4) 


0.06(1) 


2.8(2) 


9.5 


4.8 


2vp+ 




3123.30(4) 


0.06(1) 


4.8(2) 


7.7 


4.3 


V,-ot 




3.985(7) 











Notes. Phase d is given according to the definition cos {2nvt - 0) 
with respect to the epoch HJD 2454523.29297. v^at was derived from 
a weighted average of the spacings of the main septuplet. 



is Vp = 1557.67(1) [xHz, which is closer to what [Handler et al. 



found for their data set. If this value is changed to the exact 
previously known value, the residuals to the fit show no signifi- 
cant change. Therefore, we also fixed the pulsation frequency to 
Vp = 1557.6529 [,iHz. We then proceeded with a pre-whitening 
sequence to look for additional frequencies. After each iteration 
we allowed the multi-sine fit to improve the frequencies, am- 
plitudes, and phases for any additional signal found, as well as 
amplitudes and phases of all nonuplet frequencies. I.e., we kept 
Vp and Vrot fixed to the values matching the results of Handler] 
et al. The outcome of this analysis is presented in Table|4]aria 
a comparison of the forced and the unconditional solutions is 
shown Fig. [7] 
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Table 4. Multi-frequency fit results for HD 99563, forcing t he 
frequency solution to match Vrot given in Handler et al. ( 2006| l. 



id 




frequency 


amplitude 


e 


SNR 






[^Hz] 


[mmag] 


[rad] 




Vl = 


Vp 


1557.6529 


0.23(1) 


4.35(4) 


14.9 


V2 = 


Vl 


1553.678(2) 


1.62(1) 


2.619(6) 


105.1 


Vi = 


Vp+I 


1561.628(6) 


0.41(1) 


3.73(2) 


26.1 


V4 = 


V2 


1549.70(4) 


0.06(1) 


3.0(2) 


3.9 


V5 = 


Vp+2 


1565.60(3) 


0.07(1) 


0.1(1) 


4.7 


V6 = 


Vp-3 


1545.73(2) 


0.15(1) 


0.87(6) 


10.1 


Vl = 


Vp+3 


1569.58(2) 


0.10(1) 


0.6(1) 


6.1 


Vg = 


Vp-4 


1541.8(1) 


0.02(1) 


4.8(6) 


1.1 


yg = 


Vp+4 


1573.6(2) 


0.01(1) 


3.7(8) 


0.8 


VlO 




1561.671(3) 


0.94(1) 


5.32(1) 


59.5 


Vll 




1553.38(1) 


0.17(1) 


3.80(6) 


11.2 


Vi2 




1554.08(1) 


0.14(1) 


5.95(6) 


9.4 


Vi3 = 


Tcand 


1555.85(4) 


0.06(1) 


0.7(2) 


3.9 


2vp 




3115.36(2) 


0.13(1) 


2.76(8) 


9.3 


2vp- 




3107.46(3) 


0.09(1) 


0.4(1) 


6.9 


2vp - 




3119.38(4) 


0.06(1) 


2.8(2) 


4.8 


2vp+i 


3123.29(4) 


0.06(1) 


4.8(2) 


4.3 


V,ot 




3.9749 









Notes. Frequencies which are not part of the multiplet or harmon- 
ics have been labelled in order of decreasing SNR. The labels do not 
necessarily correspond to the frequencies in Table[3] Phase 6 is given 
according to the definition cos {2nvt — 0) with respect to the epoch 
HJD 2454523.29297. 



Table 5. Frequency analysis results for HD 134214. 



id 


frequency 


amplitude 


e 


sig 


SNR 




[|xHz] 


[mmag] 


[rad] 






Vl 


2949.537(1) 


1.820(8) 


0.744(4) 


4291.2 


170.1 


2vi 


5899.12(2) 


0.083(8) 


0.97(9) 


25.4 


8.2 


V2 


2915.70(1) 


0.174(8) 


5.13(4) 


106.2 


15.6 


V3 


2779.52(1) 


0.157(8) 


0.26(5) 


88.2 


13.2 


V4 


2983.33(2) 


0.116(8) 


5.81(7) 


48.4 


10.2 


V5 


2646.96(3) 


0.063(8) 


3.2(0.1) 


14.8 


5.4 


V6 


2721.96(3) 


0.061(8) 


3.2(0.1) 


14.0 


5.8 


w 


2805.31(4) 


0.054(8) 


2.0(0.1) 


10.6 


4.3 




2841.63(4) 


0.051(8) 


4.2(0.1) 


9.5 


4.3 


V9 


2687.18(4) 


0.048(8) 


0.3(0.2) 


8.5 


4.3 


VlO 


2699.84(4) 


0.044(8) 


0.8(0.2) 


7.1 


3.9 



Notes. Phase 6 is given according to the definition cos {Invt - 6) with 
respect to the epoch HJD 2454575.045398. 
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Fig. 6. Same as Fig.jijbut for HD 134214. Notice that the har- 
monic to the main frequency can be seen as a small peak at 5900 
[iHz. 
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Fig. 7. Schematic pulsation spectrum of HD 99563 (top) and 
HD9289 (bottom). Labels indicate multiplet frequencies and 
coiTesponding combinations. The thin lines show the residual 
amplitude spectra after pre-whitening of the identified frequen- 
cies. The top panel shows both the unconditional (black circles) 
and forced (grey squares) solutions for HD 99563. Note the dif- 
ferent scales for both panels. 



3.3. HD 134214 

The case of HD 134214 is more straightforward, since we can 
rule out rotational modulation on timescales of our observations. 
As shown in Fig.|6] the data reveal a main pulsation frequency 
at 2949.5 (xHz with a large amplitude of 1.8 mmag. Due to the 
strength of this signal the same instrumental modulation effect as 
in HD 99563 can be observed to even higher integer multiples. 
To be most conservative in our analysis, we therefore excluded 
any other signal that is close to vi ± «Vmost for any n. Naturally, 
we also excluded the usual integer multiples of the MOST or- 
bital frequencies as well. Applying these criteria to the whole 
range considered (down to zero and up to 6000 \iHz), we retain 
10 significant frequencies below 3000 iiHz. These are shown in 
Fig.[8] Their values are listed in Table|5] 



4. Discussion 

4. 1 . The pulsation of HD 9289 

The results presented here for HD 9289 are inconsistent with the 



findings of Kurtz et al. ( 1994j). Their analysis yielded three fre- 
quencies, vk,i = 1585.06(1) (xHz, vk,2 = 1554.79(1) (xHz, and 
i^K,3 = 1605.39(1) (xHz. The authors acknowledged, however, 
that due to severe Id"' aliasing the frequencies might have been 
confused with their aliases. A comparison to Table 2] shows that 
the primary frequency is similar in both data sets but the exact 
values do not agree within the uncertainties. We do not find any 
frequencies that have values close to vk,2 or vk,3, and we do not 
expect any instrumental artefacts that would possibly complicate 
the detection of these frequencies. However, vk,3 is consistent 
with an alias of V2, with vk3 - V2 + 2c within the uncertainties. 
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Fig. 8. Schematic pulsation spectrum of HD 134214. Squares in- 
dicate the frequencies reported in Kurtz et al. (2007). The error 
bars show our estimated upper limits for their uncertainties de- 
rived from their data set length. The arrow in the upper right 
corner indicates that the main pulsation frequency has a much 
higher amplitude (=s 2 mmag) which would explode the scal- 
ing of the figure. The residual amplitude spectrum after pre- 
whitening of the identified frequencies, and a 50-point running 
average, is also shown. 
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Fig. 9. Time-resolved analysis of HD 9289 
tiplet. The variation of the mean stellar brig 
bottom panel for comparison. 
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Here, c « 11.574 is the conversion factor from d ' to ^iHz. A 
possible match for vk,2 ~ -2c is much more strenuous and the 
disagreement is too large to reconcile it with the uncertainties. In 
any case, it is plausible that the results from |Kurtz et al.| ( |l994| l 
suffered a misidentification due to aliasing. It is not impossible, 
however, that the different results stem from an intrinsic change 
in HD 9289's pulsation behaviour 

We know from the rotational modulation of the light curve 
that Viot ~ 1.35 iiHz. Thus, the spacing in the main "triplet" (vi, 
V2, V3) appears to be twice the rotation frequency. Similarly, the 
second multiplet (V4, V5, V(,) also contains information about the 
rotation in its spacings. Whatever the nature of these modes is, 
if our interpretation is correct, this allows us to obtain an es- 
timate of the rotation frequency which is independent of the 
mean variations induced by spots. Using 1/crJ as weights, where 
(Ty are the Icr frequency uncertainties derived as described in 
Sect.[3] we calculated the weighted average of the spacings of 
adjacent peaks in the multiplets. Note that this corresponds to 
a weighting with the spectral significan ce (or the square of th e 
SNR), since crl oc 1/sig cc l/SNR^ ( [Kallinger etaL]|2008| l. 
The result is Viot.puis - 1.361(7) ^iHz, corresponding to a period 
f rot = 8.50(4) d. This estimate is more precise than the estimate 
from PDM. Therefore, invoking the 0PM, we actually observe 
components of (at least) a quintuplet with two components un- 
detected (vi -I- Vrot and vi - Vrot). Hence, under this assumption 
the mode cannot be described by a pure dipole. The question 
that remains is whether the undetected components are actually 
missing, or whether they are simply too small to be significant 
in the MOST data. 

Fig.|9] shows a time-resolved analysis of HD 9289's pulsa- 
tion. We fitted a cosine with fixed frequency vi but free am- 
plitude and phase to small subsets of the light curve after pre- 
whitening of all frequencies not related to the main triplet. A 
subset length of 1 d, analyzed in steps of 0. 1 d, was chosen as the 
best compromise between reducing the uncertainties of the fitted 



parameters and increasing the time resolution. The result is that 
the pulsation amplitude is modulated twice per rotation cycle, 
while the pulsation phase remains roughly constant throughout. 
Given the large uncertainties at epochs of low pulsation ampli- 
tude, there might be changes in phase that simply cannot be reli- 
ably detected with the current data. A strong dipole component, 
showing two amplitude maxima per rotation cycle connected to 
the two pulsation poles, would produce a phase shift of n at 
quadrature. The fact that we do not see a phase shift but a clear 
amplitude modulation, is consistent with HD 9289 being domi- 
nated by a quadrupolar {( - 2, m - Q) pulsation mode. 

The second mode can also be explained through its spacing 
as a multiplet with undetected components. Assuming the sim- 
plest solution, a quintuplet, the two "missing" components are 
the central component Vc and Vc + Vrot. Another possibility is 
a septuplet centred on V4 (or alternatively V4 + 2viot, which we 
could not detect) with four components missing. In any case, 
while still consistent with the OPM, explaining this multiplet re- 
quires a strong departure from the usual "distorted dipole" pic- 
ture. Using Vrot.puis derived above as an estimated value for v,ot, 
and calculating the position of the central (missing) frequency 
from a weighted mean based on the measured values of the de- 
tected multiplet components, the actual pulsation mode would 
have a frequency Vc = 1574.55(2) 

We repeated the time-resolved analysis by fitting a cosine 
at fixed frequency v^, after subtracting all frequencies related to 
the main pulsation multiplet. The results are presented in Fig. 10 
Due to the large uncertainties, the behaviour of the amplitude is 
not well defined, but it appears to be somewhat larger around the 
minimum in mean light. The phase changes continuously and 
linearly with time and jumps by 2n after about 1 rotation period. 
From the point of view of frequency analysis, this is due to the 



' Using Vio, instead, the central frequency of the second quintu- 
plet would be Vc = 1574.48(2) ^iHz 
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Fig. 10. Time-resolved analysis of HD9289's secondary pul- 
sation multiplet, assuming that the real mode has Vc - 
1574.55 [iHz. The variation of the mean stellar brightness is 
shown in the bottom panel for comparison. 
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Fig. 11. Time-resolved analysis of HD99563's main pulsation 
mode, choosing the frequency value displayed in Table[4] Note 
that the phase (as in cos (wf - 0)) has been transformed to (p (as 
in sin (tat + (p)) for comparison to Handler et al. ( 2006[). 



presence of V4 and the actual absence of v^. The time-resolved 
analysis essentially tries to compensate for picking the "wrong" 
frequency by shifting the phase continuously. However, if the 
physical mode frequency is Vc, the continuous change in phase 
is real. Alternatively, we could choose V4 as the physical mode 
frequency, thus requiring a septuplet with missing components. 
In this scenario, the phase does not show any jumps and remains 
rather constant while the amplitude modulation looks the same. 
We summarize that the SNR is too low for a conclusive time- 
resolved analysis of the second multiplet. 

If both modes are considered quintuplets, a principle inter- 
pretation is that both modes are dominated by ({ = 2, m = 0). 
However, in this case they should show similar relative ampli- 
tudes (or amplitude ratios) for the various components. This 
is clearly not the case. Different pulsation axes might explain 
this observation, but we refrain from speculation at this point. 
Several roAp stars recently observed with the Kepler satel- 
lite show comparable multiplet structures, and strong rotational 
modulation due to surface features (Balona et al. 201 1 , in press), 
or even give direct evidence for different pulsation axes (Kurtz, 



et al. 2011 in press). Hopefully, these more detailed observa- 



tions will help us to also improve our understanding of HD 9289. 



4.2. The pulsation of HD 99563 - a new candidate frequency 

Fig. TT] shows a time-resolved analysis of the main pulsation 
mode of HD 99563 after transforming the pulsation phase and 
defining the rotation phase so that a comparison with Fig. 6 from 
[Handler et al.| ( |2006[ l is feasible. As can be seen, the result is 
virtually indistinguishable. The difference in the two amplitude 
maxima is marginal, and the behaviour of the pulsation phase 
is exactly the same. This result does not depend on which of 
the two solutions for the exact frequency derived in Sect.|3.2| 



we chose. Also, the harmonics and linear combinations above 
3000 [xHz are the same within the error bars for both solutions. 

However, there are differences in the amplitudes of the most 
significant multiplet components between the forced and the un- 
conditional solution. This is due to the change in phase that 
comes with a change in frequency. As the relative phases of 
the components must change with a change in frequency, their 
relative amplitudes must also change in order to reproduce the 
light curve behaviour. Since the change in residuals between 
the forced and unconditional solution are negligible, there is no 
strong preference for either solution. However, if the values and 
uncertainties given in Handler et al.| (12006 ) are correct, and as- 



suming that HD 99563 has not changed and that the OPM is ap- 
plicable, the forced solution is to be preferred. Concerning the 
multiplet, the remaining difference between our results and those 
by Handler et al. ( 2006 ), in particular the relative amplitudes, are 
expected due to the different bandpasses used the in photomet- 
ric observations. Just as the amplitude of the rotational modula- 
tion varies from filter to filter, so does the amplitude modulation 
of the distorted dipole. The deviations from [Freyhammer et al. 
(j2009)l are also not surprising, since the radial velocities gen- 
erally show different relative amplitudes from photometry. One 
should be careful when using the amplitude ratios and phases 
to infer the general geometry of the oblique pulsator through, 
e.g., the technique of 'Kurtz ( 19921. In particular, one needs to 
take into account how uncertainties and instrumental properties 
affect the results. 

Concerning any additional frequencies independent from the 
main multiplet, a number of such candidate signals were for- 
mally significant in both frequency solutions. Most of the candi- 
dates in the unconditional solution cannot be exactly reproduced 
by the forced solution, but instead, other significant frequencies 
with similar but different values emerge. The only exceptions are 
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Vi2 (unconditional) and vn (forced), which match in frequency, 
amplitude, and phase. Thus, this seems to be a signal which is 
independent of the exact form of the main multiplet. We there- 
fore conclude that there is an additional candidate frequency at 
Vcand = 1555.84(4) [iHz. It is plausible that this frequency could 
not be detected by [Handler et al. since its amplitude is compa- 
rable to that of Vp_2 and Vp+2 which they were also unable to de- 
tect. We have no explanation for the remaining signals, formally 
significant but different in each solution, but we can speculate 
that they are artefacts caused by either an intrinsic or instru- 
mental long-term modulation of the main pulsation frequency. 
This would explain why their exact values are affected by the 
choice of the frequency solution, while Vcand is unaffected. The 
fact that no modulation was observed by Handler et al. increases 
the probability of an instrumental origin. Finally, in the forced 
solution, the outermost components of the nonuplet are not sig- 
nificant with a SNR of 1 . 1 and below. We can therefore not sub- 



stantiate the claim of a nonuplet by |Freyhammer et al. ( 2009, ) 



4.3. The pulsation of HD 134214 - a twin of HD 24712? 

As presented in Fig. [8] our results of a very conservative analy- 
sis for HD 134214 show that we can confirm most frequencies 
reported by [Kurtz et al.| ( 2007| l. There is a more or less regu- 
lar frequency spacing that can be found throughout the region 
between 2600 and 3000 (xHz. Only two of the frequencies we 
found do not fit this regular pattern: V3 and viq. To extract the 
most probable frequency spacing, we created a synthetic am- 
plitude spectrum in which each frequency was represented by 
a Gaussian. The width of the Gaussians was determined us- 
ing the estimated frequency uncertainties. The auto-correlation 
of this spectrum shows several peaks, but the most prominent 
peaks are a double peak at 33.80 ± 0.05 [xHz, and a single broad 
peak at 67.63 ± 0.03 \iHz. Interpreting the latter as the large 
frequency separation, we created the Echelle diagram shown 
in Fig. 12 It immediately reveals the striking regularity in the 
spacing of the frequencies. Most of the frequencies line up at 
V mod 67.63 ^iHz s; 10(,iHz, but the primary frequency vj does 
not! Considering only the frequencies that are consistent with 



Kurtz et al. (2007), vi is the only mode that clearly seems to re- 
late to a different spherical degree from the remaining frequen- 
cies. Moreover, two of the new candidate frequencies seem to 
continue the ridge of vi towards lower frequencies. This strongly 
suggests that HD 134214 pulsates in at least two different dom- 
inant spherical degrees. Kurtz et al. ( 2007| l could not yet decide 
between stable and unstable frequencies, and why most of them 
could not be observed in ground-based photometry. However, 
from the point of view of the long-timebase MOST observations, 
the frequencies seem to be stable, and HD 134214 reveals an im- 
pressive p-mode spectrum. 

To put this result into context, the value found for the large 
frequency separation is very similar to that of the well-studied 
roAp star HD 24712 (see, e.g., [Saio et al.pOlO ). Moreover, the 



frequency region of pulsation is also approximately the same. 
Taking the estimate of the stars' fundamental parameters from 
Ryabchikova et al. ( 2007 1, both are furthermore comparable 
in terms of magnetic field and effective temperature. However, 
HD24712 does not exhibit such large differences in amplitudes 
between its individual modes, shows rotational signature in its 
pulsation spectrum, and is found to have a higher v sin /. 

Putting these bits of information together allows for an in- 
teresting interpretation of the pulsation in HD 134214. At higher 
magnetic field strengths the pulsation amplitude distribution for 
different modes can become quite complex (e.g., see Fig. 12 
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Fig. 12. Echelle diagram for the frequencies extracted from the 
MOST data of HD 134214. Frequencies consistent with prior ob- 
servations of Kurtz et al.[ ( |2007| l are shown as filled black circles. 
New candidate frequencies are shown as shaded grey circles. 
Other symbols correspond to a reference model (see text) with 
1 defined as in Saio et al. ( 2010| . The frequency uncertainties 
are smaller than the symbols used. Lines are visual aids to help 
separate the two possible mode ridges. 



in [Gruberbauer et aLl[2008[ ). This can affect our ability to de- 
tect such modes. Therefore, different dominant spherical degree 
components (as can be inferred from Fig. 12 1 and a very small in- 
clination angle and/or obliquity (as can be inferred from the ab- 
sence of rotational features and the comparison with HD 24712) 
could conspire to produce the large contrast in the mode ampli- 
tudes. If the mode geometry of vi is strongly dipolar and con- 
centrates most of its pulsation power very close to the rotation 
axis, a low inclination would result in a large observed pulsation 
amplitude. On the other hand, if the pulsation of the modes on 
the left ridge in Fig. [12] is restricted to the mid-latitudes or the 
stellar equator, low inclinations will make those more difficult 
to detect. The pulsation amplitude maxima would be observed 
near the stellar limbs, and, due to their non-dipolar nature, they 
would also be subject to cancellation. This would also explain 
why [Kurtz et"aL find the pulsation amplitudes to increase to- 
wards lines which are higher up in the atmosphere, since those 
are more dominant towards the stellar limbs. 



As a reference. Fig. 12 also shows a magnetic model pro- 
duced by H. Saio (priv. comm.) using the method outlined in 
Saio et aT (2010 1. The model parameters are M = 1.65 Mg, 
logR = 0.2694, logTeff = 3.8584, logL = 0.9250, X = 0.70, 
Z = 0.02 and Bp - 4.1 kG. Although the model does not ex- 
actly correspond to our "two ridge" hypothesis, it supports the 
similarity to HD 24712. It also identifies the main mode as being 
dominated by Z = 1, while most other modes are dominated by 
higher-degree spherical harmonics. In the case of small inclina- 
tion and small obliquity, this is in support of our explanation for 
the large amplitude difference between vi and the other modes. 
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However, there are two problems with the picture presented 
above. Firstly, it cannot explain the strong amplitude contrast 
between vi and the other two frequencies on the right ridge, 
and vg. The mode spectrum of HD 24712 reveals that, over two 
radial orders, the amplitudes of modes of the same dominant de- 
gree can decrease by about a factor of 5. This is not enough to 
explain our amplitude ratios of ~ 35. Therefore, even invoking 
our hypothesis, there is still need for rather large differences in 
the intrinsic amplitudes and/or different dominant components. 
Secondly, current models suggest that for a star like HD 24712, 
the amplitude of the radiative flux variation due to pulsation has 
its maximum at the magnetic poles, regardless of the dominant 
value of / (see Fig. 12 in Saio et al. 2010). Therefore, the ob- 
served amplitudes should not have a large dependency on as- 
suming comparable intrinsic amplitudes. If both stars are indeed 
very similar, and under the assumption that the models are appli- 
cable, the scenario given above is therefore unlikely to be true. 

Whether or not our interpretation of HD 134214 is correct, 
the situation offers exciting potential to further study the simi- 
larities and differences between HD 134214 and HD 24712. As 
a first step, asteroseismic modelling of HD 134214 will hope- 
fully allow us to constrain the different mode geometries from 
a theoretical perspective. This will be presented in an upcoming 
paper (Saio et al., in preparation). 

5. Conclusions 

In this paper we have presented three data sets of roAp stars 
obtained with the MOST satellite. All three data sets have 
revealed previously unknown features, or helped to clarify 
ambiguities that could not be resolved from the ground. 

HD9289: 

• We were able to derive the star's rotation period for the first 
time withProt = 8.50(4) d. 

• At the time of observation, HD 9289 pulsated in two modes, 
and showed strong amplitude variations resulting in fre- 
quency multiplets consistent with the OPM. 

• Both modes have different multiplet structure, hinting at dif- 
ferent mode geometries. 

HD 99563: 
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The MOST observations present a picture mostly consistent 
with previous results. However, we identified a new candi- 
date frequency independent of the primary multiplet. 
We were unable to confirm the presence of a nonuplet as 
claimed by Freyhammer et al. ( 2009 1. 



HD 134214: 



MOST reveals a rich spectrum of p modes consistent with 
findings by Kurtz et al. ( 2007) l. Their stability was es- 
tablished, and the frequencies are now known to much 
higher precision. This turns HD 134214 into one of the most 
promising roAp candidates for asteroseismology. 
MOST also detected three new candidate frequencies, of 
which two fall nicely into a well-defined p-mode pattern. 
HD 134214 was found to show pulsation frequencies and 
a large frequency spacing similar to the well-studied star 
HD24712(HR1217). 
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